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Spin relaxation under identical Dresselhaus and Rashba coupling strengths in GaAs quantum 
wells is studied in both the traditional collinear statistics, where the energy spectra do not contain 
the spin-orbit coupling terms, and the helix statistics, where the spin-orbit couplings are included 
in the energy spectra. We show that there is only marginal difference between the spin relaxation 
times obtained from these two different statistics. We further show that with the cubic term of the 
Dresselhaus spin-orbit coupling included, the spin relaxation time along the (1,1,0) direction becomes 
finite, although it is still much longer than that along the other two perpendicular directions. The 
properties of the spin relaxation along this special direction under varies conditions are studied in 
detail. 

PACS numbers: 72.25.Rb,72.15.Lh,71.10.-w,67.57.Lm 
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I. INTRODUCTION 

Semiconductor spintronics which aims at making use 
of electron spin degrees of freedom has attracted much 
attention both theoretically and experimentally 1 ' 2 owing 
to its potential applications such as spin transistors. 3. 
The realization of such device requires a long spin re- 
laxation/dephasing (R/D) time which should be at least 
longer than the time of one expected operation. There- 
fore it is essential to understand the mechanism of the 
spin R/D and to find effective ways to manipulate the 
spin relaxation time (SRT) and/or spin dephasing time. 

It is well known that in n-type Zinc-blende semicon- 
ductors, the leading spin R/D mechanism is the so called 
D'yakonov-Perel' (DP) mechanism, 5 which is due to the 
spin state splitting of the conduction band at k ^ from 
the spin-orbit interaction in crystals without inversion 
center. This splitting is equivalent to an effective mag- 
netic field (EMF) h(k) acting on the spin, with its mag- 
nitude and orientation depending on k. In quantum wells 
(QW's) this EMF is composed of the Rashba term 6 due 
to the lack of the structure inversion symmetry and the 
Dresselhaus term 7 due to the lack of the bulk inversion 
symmetry. Their contributions to the energy spectra 8,9 
as well as the spin R/D have been widely studied 
in the literature. 9,10,11,12,13,14,15,16,17,18,19,20,21,22,23 Un- 
der some special conditions, such as when the growth 
direction of the QW is along (110) axis 16,24,25 or 
when the growth direction is along (001) axis but the 
strengths of the Rashba and the Dressclhauls terms 
are comparable, 4,8,9,26,27,28,29,30 the SRT shows strong 
anisotropy and along some special polarization direction 
one can get extremely long SRT (It can be infinite when 
the cubic terms in the Dresselhaus term is ignored and 
the strengths of the Rashba and the Dresselhauls terms 
are identical). Schliemann et al. further proposed a new 
spin transistor based on this strong anisotropy of the SRT 
when the strengths of the Rashba and the Dresselhauls 



terms are identical. 4 Very recently Jiang and Wu pro- 
posed to use strain to get extremely long spin dephasing 
time (also relaxation time) in GaAs (001) QW's where 
there is no anisotropy along different directions. 31 

Among these studies, Wu proposed a many-body ap- 
proach by setting up the kinetic spin Bloch equations 11 ' 12 
and solving it self-consistently with all the scattering, 
i.e., the electron-electron, electron-phonon and electron- 
nonmagnetic impurity scattering explicitly included. By 
this approach, one not only gets the spin R/D time due 
to the single-particle effective spin-flip scattering, but 
also obtains that due to the many-body effects 18,19,32 
which have shown to be dominant in some n-type semi- 
conductor QW's 18,19,33 and give many different behav- 
iors from the single-particle spin R/D both in spin 
precession 18,19,33 and spin transport. 17 The many-body 
effect here is not only referred to be the effect of the 
Coulomb scattering, but also expanded to include the 
spin R/D due to the inhomogeneous broadening from the 
EMF 18,19,32 combined with the spin conserving scatter- 
ing. Moreover, this approach also includes the counter- 
effect from the scattering to the inhomogeneous broad- 
ening from the anisotropic EMF, which has been shown 
to be also very important in the spin kinetics. 18,19 In this 
paper, we apply this approach to study the spin relax- 
ations in GaAs (001) QW's with the identical Rashba 
and Dresselhaus spin-orbit coupling strengths, especially 
at the case when the cubic term of the Dresselhaus term is 
included. The later has been investigated in the single- 
particle picture 9,28,29 without counting the many-body 
effects. 

Another issue we address in this paper is the influ- 
ence of the two spin eigenstates on the spin kinetics. It 
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as 



is noted that in our previous works 
well as many other works 9,10,13,15,20,21 in the literature, 
the equilibrium state is taken as the Fermi distribution 
of electrons in conduction band without the spin-orbit 
coupling from the DP term. Therefore, the energy spec- 
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trum of electrons is Ek — k 2 /2m* , with m* denoting the 
effective mass of electrons, and the eigenstates of spin are 
the eigenstates of a z , i.e. xt = (1) 0) T an d Xl = (0, 1) T , 
which are collinear in the laboratory coordinates. In the 
following we refer the Fermi distribution composed by 
these eigenstates (without the spin-orbit coupling) to be 
the collinear Fermi distribution and the statistics to be 
the collinear statistics. In the meantime, there is another 
equilibrium state used in the literature which is the Fermi 
distribution of electrons with the eigenstates being those 
of electrons in the conduction band with the DP spin- 
orbit coupling. Therefore the energy spectrum is now 
being 



e Kn = k 2 /2m* + V \h(k)\ 



(1) 



with rj = ±1 for the two spin branches. The eigenfunc- 
tions of spin are now 



1 . h{k) . 



(2) 



with/i(k) = h x (k)+ih y (k). The spin polarizations there- 
fore strongly depend on the momentum k. In the follow- 
ing we refer the Fermi distribution composed by these 
eigenstates (with the spin-orbit coupling) to be the helix 
Fermi distribution and the statistics to be the helix statis- 
tics as the eigenstates shows helix spin structure when the 
DP EMF is composed only by the Rashba term. 34 Some 
works are performed with the helix statistics especially 
those dealing with spin relaxation in quantum dots 35,36 
and spin Hall conductivity. 37,38 ' 39 Very recently Lechner 
and Rossler 40 and Grimaldi 41 gave some formal theory 
to discuss the spin R/D using the helix statistics. How- 
ever, whether the collinear statistics widely used in the 
spin R/D calculation is adequate in two dimensional elec- 
tron gas (2DEG) is unknown yet. In principal, when the 
DP spin-orbit coupling is weak, the collinear statistics is 
good enough for the spin dephasing calculation. How- 
ever, if one comes to the strong coupling regime which 
can be obtained in 2DEG confined in narrow quantum 
wells, 42 ' 43,44 ' 45 then two things happen: One is that the 
strong DP EMF causes strong interference effect and in- 
duces strong spin depolarization which is irrelevant to 
the spin dephasing effect; Another is that the statistical 
equilibrium state might have to be changed into that in 
the helix statistics. The later then changes all the scat- 
tering and might change the spin R/D when the DP EMF 
is strong. Nevertheless, it is not so obvious as when the 
EMF is very strong, the spin polarization decreases dra- 
matically due to the interference effect long before the 
scattering plays its role. In this paper we address this is- 
sue by solving the kinetic spin Bloch equations with both 
collinear and helix statistics. 

We organize the paper as follows: We present our 
model and kinetic spin Bloch equations with both the he- 
lix and the collinear statistics in Sec. II. Then we present 
our numerical results in Sec. III. We first compare the 
spin relaxation with two spin statistics in Sec. Ill A and 



IIIB, but without the Coulomb scattering. Then we dis- 
cuss the spin relaxation under the identical Dresselhaus 
and Rashba coupling strengths in GaAs QWs with the 
collinear statistics by including all the scattering. We 
conclude in Sec. IV. 



II. MODEL AND KINETIC SPIN BLOCH 
EQUATIONS 

We construct the kinetic spin Bloch equations 11,12 by 
using the noncquilibrium Green function method 46 with 
these two different statistics and write them in the unified 
form as follows: 



Pk + Pk|coh + Pklscatt 







(3) 



Here pk represents a single particle density matrix of elec- 
trons with wavevector k. Often we project the density 
matrix in the laboratory spin space (the collinear spin 
space) which is constructed by basis with spin states ori- 
ented relative to a fixed direction, i.e., \] an d Xl- In such 
spin space, the diagonal matrix elements Pk,cr,<7 = ,/k.cr 
with a =\ or J. describe the distribution functions of 
electrons with momentum k and spin a orientating along 
the laboratory z-axis. The off-diagonal elements pj^ de- 
scribe the corrections (coherence) between the spin-up 
and spin-down states, with its real and imaginary parts 
standing for the spin polarizations along the x- and y- 
axes respectively. 

The coherent part of the two sets of Bloch equations 
with different statistics are the same and are given in the 
collinear spin space as 



Pklcoh = i[H R + H D - ^2 VqPk- 



PkJ 



(4) 



with H R r D \ = h R r D \ ■ er representing the Rashba 6 
(Dresselhaus 7 ) Hamiltonian and erstanding for the Pauli 
matrices. It describes the spin precession along the di- 
rection of the EMF. It is noted that the Hartree-Fock 
term lq/5k-q in the coherent part cannot flip the to- 
tal electron spin as 



(5) 



Therefore the DP term is the only one that flips the total 
spin polarization. For GaAs (001) QW with well width 
a and under the infinite- well-depth assumption, h^(k) = 
a(k y ,-k x ,0) and h D (k) - 7 (fc x [fc 2 - (^) 2 ], k y [(^) 2 - 
k 2 ], 0). Here a is the strength of the Rashba term and 7 is 
the material-specific strength of the Dresselhaus term, a 
can be tuned by the external gate voltage. 6 By changing 
the external gate voltage, one may have (j = 7(f) 2 - 4 ' 30 
If one further ignores the cubic term in ho(k), at the 
special polarization direction (1,1,0), Eq. (2) becomes 
\rj) = (l,fy) T which is independent on k. Vq in Eq. (4) 

reads F q = £^ Ko( Jlf i+K ^ \I(iQz)\ 2 - Here k denotes 
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the static electric constant and k 2 = 4we 2 N e /(akBT) 
stands for the Debye-Hucke screening constant. N e 
represents the 2D electron density. The form factor 



it 4 sin 2 y/[y 2 {y 2 - tt 2 ) 2 ] with y = q z a/2. The 



square bracket [A, B] = AB — BA is the commutator. It 
is noted that pk|coh is independent on the specific spin 
statistics and can be expanded as 

/k,tlcoh = -Aalcoh = 2Im[PkPk,Ti] , (6) 
Pk.ulcoh = -*-Pk(/k,T - /k,i) - *^kPk,Tt i (?) 



with P k = (h*(k) - E q ^qPk- q , U ) and A k - 

EqKi(/k-q,t - /k-q,i)- 



The scattering terms include the contributions 
from the electron-non-magnetic-impurity scattering, the 
electron-phonon scattering and the electron-electron 
Coulomb scattering: 



Pklscatt = <^ \Uq\ 2 Y 5 ( £k -q.«l ~ £k^ 2 )T k - q , 4l (pk - pk-q)T k ,£ 2 

+ 71 Y l^^l 2 Y 7k_ q , £l {5(e k _ q;€l - £ k ,£ 2 + ^q ?z A)[(-/Vq 9a A + 1)(1 - Pk-q)pk 

q<jzA £1,62 

- A q< ^ A p k _ q (l - Pk)] 

+ ^(£k-q,6 - £k,6 - ^q^A)[A r q^A(l - Pk-q)pk - (Nqq z X + l)Pk-q(l - Pk)]}Tk,fr 
+ ^Y V l Y 5 ( £k '.«3 - £ k'-q,£ 4 + £k-q,6 - £k,£ 2 )7k_ q;?1 



+ 



qk' Ci ,$2 ,43,44 

x{Tr[T k / ii3 (l - Pk')pk'-qTk'-q,£ 4 ](l - pk-q)Pk 
~ Tr[T k / ;?3 p k /(l - Pk'-q)Tk'-q,£ 4 ]Pk-q(l - Pk)}?k,£ 2 j 



(8) 



in which {• • • }t is the Hermite conjugate of the same 
terms in the previous {}. The subscript £ denotes the 
spin branch which is a =t,| in the collinear statistics 
and 77 = ± in the helix statistics. in Eq. (8) is 
the impurity density and |J7 q | 2 = J2q z {^Z^e 2 /[Ko(q 2 + 

g 2 )]} 2 |/(«<7 z )| 2 is the impurity potential with Zi denot- 
ing the charge number of the impurity. |.9 q <j z A| 2 and 
Nqq z \ = [aqy(flq qz \/ ksT) — are the matrix element 
of the electron-phonon interaction and the Bose distri- 
bution function with phonon energy spectrum Q mz \ at 
phonon mode A and wavevector (q, q z ) respectively. It 
is noted that the scattering terms are different for differ- 
ent statistics. With the collinear statistics, the matrix 
T k ^ = Xk,<r is 1/2, a constant number, and the energy 
spectrum £k,£ = £k.r/ = £k- Then the scattering terms are 
exactly the same as those in Refs. 18,19. However, with 
the helix statistics, the scattering terms become more 
complicated. T k .,t = T^ v is a 2 x 2 matrix and becomes 
k-dependcnt: 



1 r h(k) 

2 [1 + "lhM 



(9) 



Moreover, the energy spectrum in the (5-functions in Eq. 
(8) £k,r; is anisotropic in k [Eq. (1)]. When the strength 
of the DP term tends to zero, the scattering terms in 



these two statistics become identical. 

Besides the collinear spin space, there is another spin 
space (the helix spin space) used in the literature 34 ' 40 
which is constructed by spin states |+) and |— ). The 
density matrix and the two sets of kinetic spin Bloch 
equations with two different statistics can also be written 
in the helix spin space and the physics itself is the same 
as that in the collinear spin space. They are given in 
Appendix A. 

It is noted that 



Yp k \ 

k 



= 



(10) 



for both statistics as all the scattering here is the spin- 
conserving scattering. 

The scattering terms tend to drive the out-of- 
equilibrium system back to the equilibrium one po(k) = 

{exp([ff (k) - fi]/k B T) + l}- 1 with H (k) = JjL for 

the collinear statistics and H (k) = + h(k) • er for 
the helix statistics. This can be seen from the fact that 
Pklscatt = when pk = Po(k) for each corresponding 
statistics. The equilibrium state po can be further writ- 
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ten as 



po(k) = [exp((H - /j,)/k B T) + 

= J2 - ^« 



(ii) 



for both statistics, with f(x) = [exp(x/fcsT) + the 
Fermi distribution and /z, the chemical potential. For the 
collinear statistics, Eq. (11) returns to the familiar case 
Po(k) = f(e k - m). 

It is further noted that besides the energy spectrum 
£k,£ m the two statistics is different, the spin polariza- 
tion in the equilibrium state S k = Tr[/9o(k)er]/2 for each 
k also differs in the two statistics. For the equilibrium 
states in the collinear statistics, there is not any polariza- 
tion S k = for each k. However, in the helix statistics 

S k = $k7T[/( £ k,+ -Ai)-/(£k,--M)] gives the spin polar- 
ization along the direction of the EMF of the DP term. 
But the total spin S° = £ k S k is still zero. 

Although in different statistics, the way of the spin 
relaxation might be different due to the different equi- 
librium states. Here we first show analytically that the 
conclusion that when the strengths of the Rashba and 
the Dresselhaus spin-orbit couplings are identical and the 
cubic term of the Dresselhaus term is ignored, the SRT 
along the direction hi = ^(1, 1, 0) is infinite, 4 does not 
depend on the statistics. This is because with the iden- 
tical strengths of the Rashba and the Dresselhaus spin- 
orbit couplings and in the absence of the cubic term of 
the Dresselhaus term, the EMF is now proportional to 
ni: h(k) = |h(k)|ni and hence Hr + Hu = |h(k)|fti • er. 
The total spin polarization along ni in the collinear space 
is S &1 (t) = S ■ ni, with S = J2k s k and S k = Tr[p k cr]/2. 
Then from Eq. (3) one has, in both statistics, 



d_ 

dt 



S Al +iJ2 Tt (I h R + H D ,Pk\n! ■ <r) = 



(12) 



with the help of Eqs. (5) and (10). It is easy to verify 
that the second term of the equation is zero as Tv([H R + 
i?D,Pk]ni-<r) = |h(k)|Tr[ni-erp k ni-CT-p k ni-<rni-er] 
I). Then the polarization Sfi^i) = S& 1 (t = 0) does not 
relax in both statistics. 



III. NUMERICAL RESULTS 

It is noted that when the cubic term of the Dres- 
selhaus term is included, the electron spin polarization 
along any direction relaxes, even when the strengths of 
the Rashba and Dresselhaus terms are identical. More- 
over, even when the spin lifetime along ni is infinity, 
the spin lifetimes along the two perpendicular directions 
n 2 = ^=(1,-1,0) and z = (0,0,1) are finite. It is not 
clear yet that how the different statistics change these 
SRT's. This is what we are going to address first in fol- 
lowing. 



It is seen that all the unknowns to be solved appear 
nonlinearly in the coherent and the scattering parts. 
Therefore the kinetic spin Bloch equations have to be 
solved self-consistently to obtain the temporal evolutions 
of the electron distribution functions /k l(7 (t) and the spin 
coherence Pk.uW- The SRT along any direction n is ob- 
tained by fitting the envelope of the quantity S& under a 
small initial polarization along n direction 



(13) 



By choosing spin dependent chemical potentials /i±i, one 
can construct an initial state with the given electron den- 
sity N e = J2ic Tr(/?k) at 4x 10 11 cm~ 2 and the initial spin 
polarization P„ = Sa/N e along the direction n at 5 %. 

In the numerical calculation, how to compute the scat- 
tering terms are the central problem. For the collinear 
statistics, this problem has been solved successfully with 
high accuracy and fast CPU speed. The details of the 
numerical scheme has been laid out in Ref. 19. For the 
helix statistics, however, things arc much more compli- 
cated. This is because in the helix energy spectrum, it 
is impossible to divide the momentum space into a set 
of discrete grids in the way that after performing the en- 
ergy conservation (^-functions) in the scattering terms, 
all the final states still belong to this set. Therefore we 
still divide the momentum space by the mesh as that in 
Rcf. f 9 but approximate the (5-function by the Gaussian 
distribution S(x) = lim cxp (— x 2 / a 2 ) / (s/ttu) and take 

(7— »0 

the width to be a = 0.5AZ? with AE representing the 
energy interval of the grids. This approximation requires 
much more grid points to converge the results. More- 
over, it also costs one more integration as the <5-function 
cannot be carried out analytically in the discrete space. 
Therefore the CPU time becomes much longer in the he- 
lix statistics. 

In the following, we study the spin relaxation under 
identical strengths of the Dresselhaus and the Rashba 
couplings and with the cubic term of the Dresselhaus 
term included. We first show the temporal evolution of 
the spin polarization along all directions, i.e., ni, h.2 and 
i in these two statistics. Then we investigate the SRT 
along the direction fii with the impurity scattering and 
the electron-LO phonon scattering included to compare 
the effect of different statics to the SRT. Finally we give 
the SRT with all the scattering including at different 
temperatures, electron and impurity densities and well 
widths with the collinear statistics. 

In the numerical calculation, we only need to con- 
sider electron-longitudinal optical (LO) phonon scatter- 
ing for the electron-phonon scattering as we focus on 
the high temperature regime (T > 120 K). The ma- 
trix element of the electron-LO phonon interaction is 
\9^.LO? = {c^/k/W + q 2 z )]}\m z )\ 2 with a = 



■5 



e 2 /j>/(2Q-^q)(k^ — Kq" 1 ). Kqo is the optical dielec- 
tric constant and £lqq z = £Ilo is the LO-phonon fre- 
quency. The material parameters of GaAs are listed as 
following: 47 m* = 0.067m with m the free electron 
mass, Koo = 10.8, k — 12.9 and Qlo = 35.4 meV, The 
strengths of the Dresselhaus SOC and Rashba SOC are 
7 = 25 meV A 3 and a — 7(f) 2 with a = 5 nm unless 
otherwise specified. 4 




0.15 

0.1 (c) 
0.05 

t ' 1 1 ' 

5 10 15 20 

t(ps) 

FIG. 1: The temporal evolutions of spin polarizations along 
(a): z; (b): 112; (c): fii in the two statistics with Af e =4x 10 11 
/cm 2 and iV; = 0.5iV e at T = 200 K. Solid curves: in the 
collinear statistics; Dashed ones: in the helix statistics. 



A. Temporal evolutions of spin polarization in 
collinear and helix statistics 

In Fig. 1, we plot the temporal evolutions of the spin 
polarization along the three directions: z [Fig. 1(a)], n 2 
[Fig. 1(b)] and fii [Fig. 1(c)] with the electron-LO phonon 
and the electron-impurity scattering included at temper- 
ature T = 200 K in the system with identical strengths 
of the Rashba and Dresselhaus coupling. The impurity 
density Ni is 0.5N e . The evolutions within these two 
statistics are plotted in the same figures for comparison. 

It is seen from the figure that with the cubic Dres- 
selhaus term included, the SRT along fii becomes finite 
(several hundred picoseconds), but still much longer than 
the other two perpendicular directions z and n 2 (less 
than 1 ps). This is because of the strong anisotropy 
from the EMF h(k) = fti(k)ni + /i 2 (k)n 2 with /ii(k) = 
V2 7 (k y - k x )[{lf + k -^\ and h 2 (k) = jk x k y !^. As 
the spin flip is determined by the component of the EMF 
which is perpendicular to the spin polarization, therefore 
different direction of the spin polarization experiences 
different EMF: for electron spin along the direction z, 
it feels the total EMF h(k); for the one along n 2 (i), it 
feels /ii( 2 )(k). For narrow QW's, the linear term is much 
larger than the cubic terms and it appears only in /ii(k). 
Therefore the spin along the direction fii feels only the 
cubic terms and hence experiences a much weaker spin 



relaxation; whereas the spin along the other two direc- 
tions feels the linear term and shows much shorter and 
similar SRT's. 

It is noted that the difference of the time evolutions 
of the spin polarizations between the helix and collinear 
statistics is marginal in all directions: The evolutions 
of spin polarization along fii are almost identical while 
those along fi 2 and z show marginal differences after 
the initial fast spin relaxation. This is because the 
spin-orbit couplings along the latter two directions arc 
much stronger than the one along fii. Nevertheless, it is 
shown that even for the very strong spin-orbit coupling 
strengths used in our calculation, the difference of the 
time evolution of the spin polarization between the he- 
lix and the collinear statistics is too marginal to affect 
the value of the SRT effectively. Therefore, for the spin 
polarized system, the statistics used in the scattering do 
not affect the SRT too much. We will further show this 
in the next subsection focusing on the spin relaxation 
along fii as we are interested in the longest SRT along 
this direction. 



B. Comparison of the SRT along direction (1,1,0) 
with helix and collinear statistics 



300 
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FIG. 2: Spin relaxation time vs. the temperature for different 
impurity densities N = 0.lN e , N = 0.5iV e and N = 1.0N e 
in the two statistics with the electron density N e = 4 x 10 11 
cm~ 2 . Solid curve: with the helix statistics; Dashed curve: 
with the collinear statistics. 

In this subsection we compare the SRT's along the spe- 
cial direction fii in the helix and collinear statistics with 
the electron-non-magnetic impurity and the electron-LO 
phonon scattering included. In Fig. 2 we plot the SRT 
as a function of the temperature for different impurity 
densities Ni = 0.1N e , 0.5N e and 1.0N e , with the elec- 
tron density N e = 4 x 10 11 cm~ 2 . The dots with solid 
(dashed) curves are the results with the helix (collinear) 
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statistics. It is seen from the figure that in both statistics 
the SRT's along n.! decrease with increase of the temper- 
ature or with decrease of the impurity density. These re- 
sults are consistent with our previous calculations, 18,19 ' 31 
as the inhomogeneous broadening induced by the EMF 
here depends cubically on the momentum k. We have 
shown in our previous work 18 ' 31 that in GaAs QW's when 
the cubic term of the Dresselhaus is dominant, the SRT 
decreases with the temperature. 

From the figure one further finds that the SRT along 
hi in the two statistics are almost the same. This is con- 
sistent with the result in the previous subsection. Thus 
it is suitable to study the SRT with one specific statis- 
tics, eg., in the collinear statistics as the scattering can 
be treated more precisely and efficiently. 



C. SRT along the direction (1,1,0) with the 
electron-electron scattering included in the collinear 
statistics 
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FIG. 3: Spin relaxation time vs. the electron density for 
different impurity densities iV; = (solid curves) and 0.5N e 
(dashed curves) at T = 120 K (♦), 200 K (•) and 300 K (▲) 
with all the scattering included. The well width a = 5 nm. 

In the previous two subsections, we draw the con- 
clusion that the SRT in the two statistics are almost 
the same. Although we did not include the electron- 
electron Coulomb scattering above as it is extremely dif- 
ficult to compute the coulomb scattering quantitatively 
in the helix statistics numerically to meaningfully com- 
pare the SRT from the collinear statistics for the rea- 
son given in the previous section, we still believe inclu- 
sion of the Coulomb scattering would not bring a vast 
change to the SRT from the two statistics, especially for 
the polarization along fii where the spin-orbit coupling 
is weak. Moreover, recent studies of spin kinetics have 
shown that the Coulomb scattering plays a crucial role 
in spin R/D. 11 ' 12,14 ' 18,19 Therefore in this subsection we 



investigate the properties of the SRT along fii in the 
collinear statistics with the Coulomb scattering explic- 
itly included. 

In Fig. 3 we plot the SRT versus the electron density 
with and without impurities at different temperatures 
T=120 K (♦), 200 K(.) and 300 K (a). The solid curves 
are for case without impurity and the dashed ones are for 
the case with impurity (TV, = 0.5N e ). We find the im- 
purity and the temperature dependence of the SRT with 
the electron-electron Coulomb scattering included is all 
the same with that in Sec. IIIB without the Coulomb 
scattering. By comparing the values at the same condi- 
tions in Fig. 2, it is easy to find that the SRT with the 
electron-electron scattering is larger, which is consistent 
with the previous results. 17,18 ' 19 . It is also seen from the 
figure that the SRT decreases with the increase of the 
electron density and the rate of decrease is slower with 
increase of the temperature. It can be understood from 
the view of the inhomogeneous broadening: 32 For higher 
electron density, electrons tend to populate to higher mo- 
mentum space. Therefore the inhomogeneous broaden- 
ing induced by the EMF becomes larger which results 
in a shorter SRT. Moreover, when the temperature is 
high enough, electrons prefer to distribute in high mo- 
mentum space even at low electron densities. Therefore 
the inhomogeneous broadening becomes less sensitive to 
the change of the electron density. This explains why 
the SRT becomes insensitive to the electron density at at 
high temperatures. 




FIG. 4: Spin relaxation time vs. 
electron densities Af e = lx 10 11 c 



the 

„-2 



width of the QW for 
(solid curves) and 4 x 



10 cm (dashed curves) at different temperature 120 K 



(curves with •) and 300 K (curves with A), 
density iV; = 0. 



The impurity 



Now we turn to the well width dependence of the SRT. 
It is noted that when the well width is changed, the re- 
lation (3 = 7(f) 2 is maintained by changing the external 
gate voltage. In Fig. 4 we plot the SRT versus the width 
of the QW for different electron densities N e — 1 x 10 11 
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cm" 2 and 4 x 10 11 cm" 2 at T = 120 K and 300 K. The 
impurity is taken to be zero in the calculation. Contrary 
to the previous results in (001) QW's with only the Dres- 
selhaus spin-orbit coupling, 18 ' 19 here one finds that the 
SRT decreases with the increase of the well width. This 
is because in the present case with identical strengths 
of the Rashba and Dresselhaus spin-orbit coupling, the 
component of the EMF which can flip electron spin along 
rii-direction is only the cubic term of the Dresselhaus 
coupling, which is independent on the well width. Nev- 
ertheless, the linear EMF along the direction hi is pro- 
portional to a~ 2 which becomes larger with decrease of 
the well width. For large linear EMF along ni, elec- 
tron spins prefer to line along ni. Therefore spin flip- 
ping along fii is suppressed and thus the inhomogeneous 
broadening which leads to the spin relaxation along fti 
is suppressed. Consequently the SRT increases with de- 
crease of the well width. 



IV. CONCLUSIONS 



Dresselhaus spin-orbit coupling. However, the SRT de- 
creases with increase of the well width, which is contrary 
to the previous results. All these behaviors are explained 
following our previous theories. 
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In conclusion, we study the spin relaxation under the 
identical Rashba and Dresselhaus coupling strengths in 
(001) GaAs QW by the kinetic spin Bloch equations with 
different statistics: the collinear and the helix statistics. 
We show that the SRT's calculated from both the helix 
statistics and collinear statistics differ marginally in all 
directions. This is very important for the study of the 
spin R/D as most spin R/D calculations are performed 
with the collinear statistics and the numerical evaluation 
of the scattering within the collinear statistics is much 
more accurate and faster. 

We show that when the cubic term of the Dresselhaus 
spin-orbit coupling is ignored, the SRT's along (1,1,0) 
direction in both statistics are infinite. However, with 
the inclusion of the cubic term, the spin relaxation along 
this direction becomes finite, although still much longer 
than that along the other two perpendicular directions. 
We study impurity density, temperature, electron den- 
sity and the well width dependence of the SRT along 
this direction. We find the SRT increases with increase 
of the impurity density and decreases with increase of 
the temperature or the electron density. These proper- 
ties are all the same with case of (001) QW with only the 



APPENDIX A: THE KINETIC SPIN BLOCH 
EQUATIONS IN THE HELIX SPIN SPACE 

In the helix spin space, the physics meaning of each 
matrix elements , = {t]\pu\v') i s n °t well defined as 
the helix spin state is a mixture of spin-up and -down 
states and is k dependent [Eq. (2)]. Here the superscript 
h denote the quantity in the helix spin space. The density 
matrices and the kinetic spin Bloch equations in the helix 
spin spaces can be transformed from that in the collinear 
spin space by a unitary transformation: p£ = U^p^Uk, 
with 

1 ( I I \ 

U± = — \ _^k)_ fe(k) • (Al) 
V2 \ | h(k) | | h(k) | J 

The coherent terms are 

/kVlcoh = -/£,_U = 2Im[i? h ^ i+ _] , (A2) 
f>i+-U = -ii£(/£ 1+ - /£,_) ^KP^ + - , (A3) 

in which 



P£ = - E W^-q] + *R^'(k, q)]Im[ Pk _ qi+ _] + ilm[fc'(k, q)](/£_ q , + - /£_ q ,_)/2} , (A4) 
q 

A£ = 2|h(k)|-£u q {Re[/ l '(k,q)](/^ (A5) 
q 
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In these equations h'(k,q) = h*(k)h(k — q)/[\h(k)\\h(k — q)|]. The scattering terms are 

PM-catt = {nNi ]T W E 5 ( £ k-q,6 - ek,fo)5 k , k _ q T^_ qi£l (5 k _ q , k pJ; - pt q S k - q , k )T k ^ 2 

+ n X] l5qg,A| 2 X 5 , k,k-q7 k _ q€l 

q<^A £i,f 2 

X {<5(£k- q , ?1 - £k,? 2 + ^qq,A)[(^ qfc A + 1)(1 - P k - q )^k-q, k Pk 

- A r qg^AP k _ q S' k _ qik (l - Pk)] 

+ £( £ k-q,£i - £ k ,C 2 - ^q 9z A)[^q 9;!i A(l - P k - q )5'k-q,kP k 

- (A^A + l)p k _ q 5 k _ q , k (l - Pk)]}T^ 2 

+ TT^Vq 2 X ^( £ k',J 3 ~ e k'- q ,£ 4 + £ k- q ,£i - £k,$ 2 )'S'k,k-q7k-q,| 1 

qk' 5l,?2^3^4 

x{Tr[T k \ ?3 (l - p^)^, k ,_ q p^_ q T k fe ,_ q;?4 5 k ,_ q , k ,](l - ^_ q )S k _ q , k ^ 
-TV[T k \ ?3 p£,S k ,, k ,_ q (l - /£_ q )T&_ qi£ >,_ q , k ,]^_ q S k _ q , k (l - ^)}T k ^ 2 | 

+ {•••} t (A6) 



with 5 k , k _ q = t/ k ?7 k _ q = ^[(1 + h'(k, q)) + (1 - space, T k£ = T kcr = 1/2 for the collinear statistics and 
/i'(k,q))a x ] and T£ >£ = U^T^U^. In the helix spin T k ^ = T k l I( = ±[1 + 77<7 Z ] for the helix statistics. 
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